Darboux transformations for quasi-exactly solvable Hamiltonians 
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Abstract 

We construct new quasi-exactly solvable one-dimensional potentials through Dar- 
boux transformations. Three directions are investigated: Reducible and two types of 
irreducible second-order transformations. The irreducible transformations of the first 
type give singular intermediate potentials and the ones of the second type give complex- 
valued intermediate potentials while final potentials are meaningful in all cases. These 
developments are illustrated on the so-called radial sextic oscillator. 



1 Introduction 

The resolution of the one-dimensional Schrodinger equation has attracted much attention 
since the first statements of quantum mechanics. It is now well known that this equation 
can be analytically solved for a small number of interactions only. In that case, it is referred 
to as an exactly solvable (ES) equation. Because of their major interest (analytic solu- 
tions describe more finely the physical reality), several methods have been proposed in order 
to increase the number of these ES Schrodinger equations. One can cite the factorization 
method introduced by Schrodinger [Ij] and actualized by Witten in the form of supersym- 
metric quantum mechanics 0. It is based on the fact that a one- dimensional Schrodinger 
Hamiltonian 

Ho - ~ + V (x) (1) 
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can be factorized as 

H = LjLo + a , (2) 

with 

L = ^ + W (x), (3) 

the constant a being called the factorization constant. One can then obtain in a straight- 
forward way the eigenfunctions of 

H 1 = L Ll + a = --^ + V 1 (x), (4) 

V 1 (x) = Vo(x) + 2W^(x) (5) 
where the prime denotes the derivative with respect to x. They are precisely given by 

<fr N (x) = L ip N (x) , (6) 

where iPn(x) stands for the eigenfunction (of eigenvalue E N ) related to H . The relation 
(f|) means that the operator Lq realizes the transformation in the opposite direction. The 
function 

ip N {x) = L\(j) N {x) , (7) 

is an eigenfunction of H differing from ^ by a normalization constant only. Note that 
neither Lq nor L\ are unitary and they do not preserve the norms of the functions. The 
relations (P) and (|7|) imply that Hq and Hi are isospectral up to the fact that either the 
ground state of H Q could belong to the kernel of the operator L or the ground state of Hi 
could belong to the kernel of the operator L* . In this case the spectra of H and Hi differ 
by the ground state level only and they are strictly isospectral otherwise. 

Another method for giving rise to ES Hamiltonians is the so-called Darboux transfor- 
mation ||. This method is based on a general notion of transformation operators J3j. The 
main element in definition of such an operator is the so-called intertwining relation 

LH = HiL . (8) 

The factorization method {a fortiori the supersymmetric quantum mechanics) appeared long 
after the publication of the Darboux paper which was unfamiliar to physicists. But since 
the paper by Andrianov et. al. [H it was realized that the methods are equivalent. The 
relation (H) gives the same expression for the transformation operator as in the original 
paper by Darboux 0] which coincides with (0)-(|]) and precises the function Wq(x) to be the 
logarithmic derivative 

= _ d{ln^x)) 
ax 

of a solution to the initial Schrodinger equation 

(H - a)ij(x) = . (10) 

Note that a here is exactly the same factorization constant as in Eqs. (|2]) and (|]). The 
function ip(x) is called the transformation function (sometimes the factorization function). 
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The only condition imposed on %p(x) is its absence of zeros inside the interval where the 
equation (pf) is solved. This condition may be satisfied only if a < E (see e.g. 0) where 
E is the ground state energy of H if it has a discrete spectrum or the lower bound of the 
continuous spectrum otherwise. For a = E the function ip(x) is evidently the ground state 
function of H and for a < E it is an unphysical solution of (|10|). 

It is evident that the Hamiltonian H\ can be taken as the initial Hamiltonian for the next 
transformation step etc. In such a way one gets a chain of so called reducible transformations. 
The main feature of such a chain is that every element of the chain has a well-defined (this 
means essentially self-adjoint) Hamiltonian acting in the same Hilbert space and the resulting 
action of the whole chain may be obtained with the help of a higher order transformation 
operator only acting on solutions of the initial equation. A natural question then arises. Is it 
possible to get something new with respect to a chain of transformations by defining an nth 
order transformation operator as an nth order differential operator satisfying the intertwining 
relation (||[)? The answer is positive: there are the irreducible chains of transformations. The 
main feature of such chains is that despite of the fact that they can also decompose into a 
superposition of first order operators, some of the corresponding intermediate Hamiltonians 
can not be defined properly: They are either singular or non Hermitian. More precisely, it 
was first established that some intermediate potentials may have poles inside the interval 
where the initial Schrodinger equation is solved but the final potential is regular || (the 
super symmetric interpretation can be found in || and a generalization of Krein's statement 
|T0||). In this paper we refer to this case as irreducible transformations of the first type. 
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Another possibility to get an irreducible chain consists in choosing complex factorization 
energies ci\ = a and a<i = a (the bar over a symbol means the complex conjugation). In this 
case some intermediate potentials are complex but the final potential is real [|1TJ. We refer 
to this case as irreducible transformations of the second type. 

Irreducible Darboux transformations have already been exploited in order to give rise 
to a large amount of new ES Schrodinger equations 0, P|-|12|. However a new kind of 



Schrodinger equations appeared in the literature some twenty years ago ]n|. These are the 
so-called quasi-exactly solvable (QES) Schrodinger equations, i.e., those for which a finite 
number of solutions can be analytically determined. These equations are rather exceptional. 
A previous attempt to list them in an exhaustive way has been performed by Turbiner |14| 
on the basis of the link between these equations and the sl(2, R) vector fields preserving a 



finite-dimensional space. However it is now understood that some of these equations [[15] do 
not have any connection with the Lie algebra sl(2, R). Therefore it is interesting to find other 
approaches giving rise to new QES equations. Supersymmetric quantum mechanics is one 
of these approaches but it has been mainly used in the context of ES equations (for a recent 
review, see e.g. |H|]). In what concerns QES equations there are only several papers on this 



subject (see e.g. |17| ). The purpose of this paper is to fill in this gap. Thus, in Sections 



2 and 3, we focus on the reducible and irreducible second-order Darboux transformations, 
respectively. In each of these Sections we produce new QES potentials starting from the 
same QES potential. For simplicity we choose it as the one corresponding to the radial 



sextic oscillator |lq even if our developments can evidently be applied to any of the already 



known QES interactions. 
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2 Reducible second-order Darboux transformations 

Let us start with the one- dimensional Schrodinger Hamiltonian Hq as given in Eq. (|]) with 

V (x) = a 2 x 6 - 2a(2M + 2s + l)x 2 + 4 ( s ~ VfKg ~ 3 / 4 ) x £ R + (n) 

ar 

which describes the radial sextic oscillator, a prototype of QES systems. In Eq. ([□]), a 
and s are positive constants while the positive integer M is related to the number (equal to 



M + 1) of analytic eigenfunctions. These functions are expressed as |T8 



i) N {x) = exp (-^V^Pjf \x 2 ) , N = 0, 1, . . . , M (12) 

where 

M 

^fV) = £cW*- (13) 

ra=0 

The coefficients have to be determined through the resolution of a system of algebraic 
equations. Precisely these coefficients are 

/ i \M—n 

C (N) - lz£i y (N) (a) (U) 



= 1 if n = M 



(M-n)\. 



V a™. if n^M, 



with the 's being known through [18 
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s-^ (4ax 4 - 8s) s-^ 4x 2 , . , 

E 1 ! — ^r- E 7i — — (aot - 4aMx - En = °- ( 15 ) 

fc=i ^ fc ^ =1 (x 2 -c^ ; )(a^-a, v ; ) 

Let us now turn to the Darboux transformation. In the reducible case we are concerned 
with, it is a product of two Darboux transformations as given in Eq. (0) with Wq(x) being 
fixed according to Eq. @. We first choose ip(x) to be the ground state of H in order to 
prove that even in this simplified context we can generate new QES potentials. Indeed with 
such a choice we obtain 

, , , (2s -1/2) 2x 
Wo(x) = ax 3 - 1 - 1 ] - ~, _(o h ■ ^ 



k=l 



\x* — a 



Therefore one can immediately obtain M eigenfunctions 4>n(x) (as expressed in Eq. (||)) of 
the Hamiltonian Hi defined in Eqs. (|j) and (Q) 



-vi = 6XPl !^f" w(PS\x 2 ),Plf\x 2 )), V 1.2 U. (17) 
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The symbol W stands for the usual Wronskian i.e. 

W(fx(x), h{x)) ee Hx)^ - h{x)^- . (18) 

Taking fi(x) and ^(z) being respectively the polynomials P^(x 2 ) and P^(x 2 ) this Wron- 
skian will be denoted by Wkn(x). It is also clear from (|T7| ) that 4>o(x) = in accordance 
with Eq. (D. 

The potential V\{x) related to these M nontrivial solutions 4>n{x) is expressed through 
Eq. © as 

-r r / \ 26 „ 1N 2 (4s 2 -1/4) ^ (2s + M-a(4 0) ) 2 ) 
^i(x) = aV-4a(M + s-l)x 2 + ^ ^ 2 7 + 4 ( V Q) k ' 1 



2 (°) 



M (0) M (0) . (0) 

'fe(* 2 -4 0) ) 2 fc fc(- 2 -4 0) )(- 2 -«! 0) )' 1 J 



It is thus a QES potential and being not subtended by sl(2, R), it has not been listed in [[14 



Yet at this stage we have thus produced a new QES potential. This is a generalization of a 
previously obtained potential |17|| which corresponds to the choice a = 1/2 and s = 1/4. 

We now proceed further by going to the second first-order Darboux operator L\ such 
that 2 

H x = -jL + Vl (x) = L\L X + Ex , (20) 

with 

Ll = ± + Wx{x) = ±- d ^ x)) . (21) 
We will then obtain once again a new QES potential V 2 (x) defined according to Eq.fl^) 

V 2 {x) = V l {x) + 2W[{x) (22) 
and the corresponding (M — 1) eigenf unctions 

X n(x) = L 1( j> N (x) , N = 2,...,M. (23) 

The function Wi(x) characterizing this second Darboux operator is determined according to 
Eqs. ( PH ) and ( ]T7| ) and it is given by 

W 1 (x) = -W (x) + (Ex - E ) M % } y\ [ ] - (24) 

It is then straightforward to obtain V 2 (x) and Xn( x ) through Eqs. ( fZ2"| ) and fl2"3"|), respectively. 
They are 

w=li (,) + ^-^ f '^y'» ) (25) 
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and 



Qj a . X 



2s-l/2 



Xn(x) = ^(-^x 4 )^ (26) 

x [-E P$(a?)W 1N (x) + E 1 P$(x 2 )W QN (x) - E N P^\x 2 )W Q1 (x)} , 
N = 2,...,M. 

Once again it is immediate to observe from Eq. (p7|) that Xo(^) — = 0- We thus 

obtain in Eq. (|2T|) the (M — 1) eigenfunctions of the new QES potential V 2 (x) given in Eq. 



As stated in the Introduction, this succession of two first-order Darboux operators is 
equivalent to a unique second-order Darboux operator defined by 

L = L X L . (27) 

Following the general theory of Darboux transformations ||, and this can be directly checked 
through Eqs. (^), fllEp, fl2~T| ) and (f?|), this operator is such that 

Lij; N (x) = W-'iM^iM^WiM^iM^^Nix)) (28) 

giving again the eigenfunctions Xn( x ) defined in Eq. (j2"T|), while 

V 2 (x) = V (x) - 2^- 2 [\nW(Mx),Mx))}- (29) 

This Darboux operator L is thus the one connecting the two QES potentials Vq(x) and V 2 (x) 
given in fl25|). The intermediate step V±(x) (see formula (0)) is also physically meaningful. 
It is precisely the statement of reducible Darboux transformations but adapted here to QES 
equations. 

We observe here the usual situation inherent for the potentials with the centrifugal term 
of the type 1(1 + l)/x 2 , (/ = 2s — 3/2). The Darboux transformation results in the change 
I — > I + 1 for this term even if it is not present in the initial potential (i.e. 1 = 0). This is due 
to zero boundary condition for the functions (O) and their specific asymptotic behaviour 
when x — > 0, iPn(x) ~ x l+1 . 

We end this Section by a specific example. This will be in particular useful in order to 
compare the results we have obtained here with the ones of the irreducible context in the 
next Section. We thus take 

a = - , M = 2 (30) 



so that the starting potential Vq(x) is 

V;M4.M5 + 2 S ). V (s - 1/4) 2 (s - 3/4) . (31) 

4 X z 
The corresponding eigenfunctions and eigenvalues are given by 

ij> (x) = exp (-\x 4 )x 2s - l/2 (x A + 2V '4s + lx 2 + 4s) , E = -A^As + 1 , (32) 
8 

Mx)= exp(-^xV s - 1/2 (x 4 -4s-2), E x = , (33) 

8 

ip 2 {x) = exp (--x A )x 2s ~ 1/2 (x A - 2V4s + lx 2 + 4s) , E 2 = 4 V / 47TT (34) 



as seen from Eq. (|j"2"|) . The new QES potential V2(x) resulting from the second-order 
reducible operator L given by (|27j) is (see Eq. ( j25[ ) ) 



V 2 (x) = -x 6 - (2s - l)x 2 + — + 



4 x 2 (x 4 + 2 v / 4sTlx 2 + 4s + 2) 

32a; 2 



(35) 



(x 4 + 2 v / 4sTlx 2 + 4s + 2) 2 

while its unique analytic eigenfunction % 2 (x) as well as the corresponding energy E 2 are 

1 x 2s + 3 / 2 
X2 (x)=exp(--x 4 ) £ 2 = 4v / 4sTT . (36) 

o (x 4 + 2\/4s + lx 2 + 4s + 2) 

Moreover there exists in this case a meaningful intermediate step characterized by the po- 
tential Vi(x), i.e., 

t^/ n 1 e 0/ iN 2 (4s 2 -1/4) 4(l + v / 47Tl) 
VAx) = -x 6 - 2(s + l)x 2 + ^ + 



4 x 2 x 2 + V4s + 1-1 

4(1 - y/isTT) 8(1 - v / 4sTT) 



x 2 + V4s + 1 + 1 (x 2 + V4S + 1 - l) 2 
8(l + v / 4sTl) 8V/47TT 



(x 2 + V4s + 1 + l) 2 x 4 + 2V4s + lx 2 + 4s 



(37) 



Two eigenf unctions and eigenvalue related to this potential are also known. They are given 
by 

j. , \ 1 1 4n 2 ,+i/2 a;4 + V4s + lx 2 + 4s + 2 

0i (x) = exp( — x 4 )x 2s+i/i — - — ; , ^i = 0, (38) 

Vy 1 PV 8 ; x 4 + 2v^Tlx 2 + 4s ' y 1 

<hix)= exp (-ix 4 )x 2s+1 / 2 % ~ 4S , J B 2 = 4 v / 4sTT. (39) 

^ V ; PV 8 ; x 4 + 2V47T1X 2 + 4s ' 1 ; 

3 Irreducible second-order Darboux 
transformations 

3.1 Transformations of the first type 

We now turn to the irreducible context. We still consider at the start the potential Vq(x) given 
in Eq. ( [TT]) as well as the corresponding eigenf unctions i[)n(x), N = 0, 1, . . . , M specified in 



Eq. (|12D . The main purpose here is to find an operator similar to the one given in Eq. (|28| ) 
but being such that the factorization Q2"7| ) cannot give rise to a physically meaningful Vi(x). 
We thus propose to consider an operator L being built on ipi(x) and ^2 0*0 i- e - 

L^(x) = ^- 1 (^(x),^ 2 (x))iy(^ 1 (x), i> 2 {x)^ N {x)) (40) 

such that LiPn{x) will stand for the (M — 1) (N = 0, 3, 4, ... , M) eigenfunctions of the new 
QES potential replacing V" 2 (x) given in (p9|). 
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One can see, comparing Eq. ( f40|) with Eq. (p8|), that the irreducible case can be deduced 
from the reducible one by the simple changes of indices: — > 1 and 1 — > 2, in Eqs. (p8|) and 
(p9|) respectively. Through example (^) we thus get 

V 2 (x) = -x b - (2s - l)x 2 + — + 



4 x 2 (x 4 - 2V4s + lx 2 + 4s + 2) 

32x 2 



(x 4 - 2V4s + lx 2 + 4s + 2) 2 



(41) 



i 2s+3/2 

X M^xp(--^) (x4 _ 2V _ x2 + 4s + 2) , Eo = -4v^TT. (42) 



and 



One immediately notice that in this irreducible case, the potential (|4"ID as well as its eigen- 
function (]42| ) can be deduced from their reducible analogues ( |35f) and (|3~6"D by a simple change 
of sign in front of the square roots. However it is a new QES potential with a unique eigen- 
function Xo{ x ) given in Eq. ( f42"]) and characterized this time by the lowest energy of the 
potential ( |TTD when restricted to the choices (]3"0"|). 

As stated before, the Darboux operator L cannot give rise to a physically relevant in- 
termediate potential Vi(x). Indeed, the operator L defined in Eq. (f40D can be written as 
L = LxLq with 

Lj = ^ + wj(x), j = 0,1. (43) 

More precisely, the function wo(x), the irreducible analogue of Wq{x) given in Eq. ([16]), is 
determined as 

, . 1 o (2s - 1/2) / 1 1 \ 

W X = -X 3 - - — - 2X ; + ; , 44 

which effectively coincides with the one of Eq. ( [IB] ) up to the choices ([30]) and the replacement 
of (k — 1, . . . , M) by af\ As clear from Eq. (|4j), this leads to singularities which were 
not present in the reducible case (|i~6|). These singularities are also recovered at the step of 
the potential Vi(x) obtained through Eq. (|5|). Indeed we have (compare with Eq. (|37])) 

-w T / \ 1 6 n , iX 2 (4s 2 -l/4) 4(1 + V¥+1) 
Vi(x) = -x 6 -2(s + l)x 2 + ^ + 



4 " x 2 x 2 - a/4s + 1 + 1 

4(1 - VisTT) 8(l + v / 4^TI) 



x 2 - v/is + l - 1 (x 2 - VAi+l - l) 2 
747+1) | Sy^sTl 



(x 2 - V4s + 1 + l) 2 x 4 - 2V4s + lx 2 + 4s 

This potential is ill-defined in the space M.q we are working in and cannot be considered as 
a convenient intermediate step proving then the irreducible character of L. 
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3.2 Transformations of the second type 

Let us now conclude with the third specific feature of the Darboux transformation namely 
the fact that one can construct physically significant new potentials from complex solutions 
of the initial equation. Once again we concentrate on the radial sextic oscillator characterized 
by the potential Vq(x) given in Eq. ( |LT|) . In order to illustrate our approach, we only need 
to find non-physical eigenfunctions corresponding to complex eigenvalues. 

We shall now show that the potential ([11]) with the integer and half-integer values of 
s has analytic solutions for some complex values of the "energy" (for simplicity we shall 
continue to use this terminology). For this purpose we note that this equation is covariant 
under the transformation s — > 1 — s, M — > M + 2s — 1. This means that for a given value 
of the positive integer M and some values of the^integer or half-integer s, the eigenfunction 
of Hq with s replaced by s = 1 — s and M by M = M + 2s — 1 is also the eigenfunction 
of Hq with the given M and s. Note that s is negative for s > 1. The values of the energy 



calculated according to Jl8| with negative s might be either complex or real. When they 
are real they lead to known solutions but for complex values of the energy we obtain new 
solutions of the Schrodinger equation. These new solutions when used to construct second 
order transformations lead to new real QES potentials. The second-order Darboux operator 
A is defined in complete analogy with Eq. (|28|) by 

Aip ( x ) = W-\^{x), ${x))W{${x), $(x), ip (x)) . (46) 

The state ([46]) will then be the eigenfunction (of the same eigenvalue as the one of ipo(x) ) 
of the QES potential 

v 2 [x) = V {x) - 2^-{hxW{^{x)^{x))) . (47) 
ax z 

Let us illustrate now this possibility by a specific example. We take a — 1/2, M = and 
s — 2. In this case M = 3 and s~ — — 1. The complex eigenfunction of the potential 

1 ^ 
V (x) = -x 6 -5x 2 + ^ (48) 
4 4ar 

is 

^( x ) = x -5/ 2 eX p(-x 4 /8)[ix 2 (ix 4 - 5) + iVh{]^x A - 1)] . 

It corresponds to the energy E = The single analytic eigenfunction (corresponding 

to E = 0) of this potential (|8|) is 

1 7 

ipo( x ) — ex P ( — x 4 )x2 . (49) 



The new QES potential reads 

4 X +X + 4x 2 ' "20 + 4x 4 + x 8 " " J (20 + 4:X 4 + x 8 ) 2 



1 3 —fir 2 -I- t 6 t 6 

= ^ + - 2 + A + 16 nn 'I If s - 2048 - ■ * , li^ (50) 



and its solution is 



A / n exp(-|a; 4 )a;2( a ; 4 + 6) 
x s + 4x 4 + 20 
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